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The current status of origin of Ultra High Energy Cosmic Rays (UHECR) is reviewed, with emphasis given to 
elementary particle solutions to UHECR problem, namely to Topological Defects and Super-Heavy Dark Matter 
(SHDM) particles. The relic superheavy particles are very efficiently produced at inflation. Being protected 
by gauge discrete symmetries, they can be long lived. They are clustering in the Galactic halo, producing 
thus UHECR without Greisen-Zatsepin-Kuzmin cutoff. Topological Defects can naturally produce particles with 
energies as observed and much higher, but in most cases fail to produce the observed fluxes. Cosmic necklaces, 
monopoles connected by strings and vortons are identified as most plausible sources. The latter two of them are 
also clustering in the halo and their observational predictions are identical to those of SHDM particles. 



1. Introduction 

Ultra High Energy Cosmic Rays (UHECR) is 
a puzzle of modern physics. Its solution needs 
the new ideas in astrophysics or in elementary 
particle physics. 
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Figure 1. AG AS A spectrum compared with "as- 
trophysical" spectrum calculated under assump- 
tions that the sources distributed uniformly in the 
Universe and have generation spectrum ~ _E~ 2 - 3 
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The problem of UHECR is known for more 
than 30 years. It consits in observation of pri- 



mary particles with energies up to 2 — 3 ■ 10 20 eV 
If these particles are extragalactic protons 
and their sources are distributed uniformly in the 
Universe, their spectrum must expose steepening, 
which starts at energy Ebb ~ 3 ■ 10 19 eV due to 
interaction with microwave photons. This steep- 
ening is known as the Greisen-Zatsepin-Kuzmin 
(GZK) cutoff §. The GZK cutoff is not seen 
in the observed spectrum. The spectrum of 
UHECR according to AGASA observations |j| is 
shown in Fig.l together with the spectrum cal- 
culated for uniform distribution of the sources in 
the Universe under assumption that generation 
spectrum is proportional to E~ 23 . The excess 
of the observed events above the GZK cutoff is 
clearly seen. The observational data for UHECR 
(E < 1 • 10 19 eV) can be summarized as follows. 

• At E > 10 19 eV the spectrum is flatter than 
at lower energies and it extends up to 2 — 
3 • 10 20 eV (maximum observed energies). 

• Chemical composition is favoured by pro- 
tons, though UHE photons are not excluded 
as primaries. 

• Data are consistent with isotropy, but close 
angular pairs (doublets) and triplets com- 
pose about 20% of all events at E > 4 • 
10 19 eV (22 events in doublets and triplets 
from 92 total j|. 
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Galactic origin of UHECR due to acceleration by 
sources located in the Galactic disc is excluded. 
Numerical simulations of propagation of UHECR 
in magnetic fields of disc and halo of the Galaxy 
predict the strong anisotropy for particles with 
rigidity E/Z > 1 • 10 19 eV (@ - |l(J). 

Extragalactic protons, if their sources are dis- 
tributed uniformly in the universe, should have 
GZK cutoff due to pion production on microwave 
radiation (see Fig.l). 

Extragalactic nuclei exhibit the cutoff at the same 
energy E <~ 3 • fO 19 eV, mainly due to e + e~-pair 
production on microwave radiation [jllj , [pj , |T^] . 

Nearby sources must form a compact group with 
large overdensity of the sources to avoid GZK cut- 
off ||. Local Supercluster (LS) with the typi- 
cal size Rls ~ 10 Mpc is a natural candidate 
for such group. The calculations (see ^) show 
that for absence of GZK cutoff the LS overden- 
sity 5ls > 10 is needed, while the observed one is 
Sls ~ 1-4 ]l3|. Note, that diffusion propagation 
due to magnetic field cannot help in softening of 
GZK cutoff. 

Nearby single source can provide the absence of 
GZK cutoff. The idea is that powerful sources 
of UHECR in the Universe are very rare and by 
chance we live nearby one of them. Such case 
has been numerically studied for the burst gener- 
ation of UHECR and their non-stationary diffuse 
propagation fl^j . Anisotropy can be small even 
at energies exceeding 1 ■ 10 20 eV. The calculated 
cutoff at E c w 1 • 10 20 eV is questioned by exis- 
tence of two events with energy 2 ■ I 20 eV. 

An interesting case of single source UHECR 
origin was recently proposed in | fl5| . The physical 
essence of this model can be explained in the fol- 
lowing way. A nearby single source is the powerful 
radio galaxy M87 in Virgo cluster. UHE parti- 
cles from it falls to gigantic magnetic halo of our 
Galaxy (with height about f.5 Mpc), where the 
azimuthal magnetic field diminishes as 1/r. Mag- 
netic field in the halo focuses the highest energy 
particles to the Sun in such way, that arriving 
particles have isotropical distribution. Numeri- 
cal simulations of the trajectories in the magnetic 
field, similar to that in the galactic wind, confirm 
this model. This interesting proposal should be 



further studied taking into account such phenom- 
ena as diffuse radio, X-ray and gamma radiation 
produced by high energy electrons diffusing from 
the Galactic disc. The calculations of these pro- 
cesses limit the size of magnetic halo by 3 — 5 kpc 

Acceleration of UHECR is a problem for as- 
trophysical scenarios. Shock acceleration and 
unipolar induction are the "standard" accelera- 
tion mechanisms to UHE, considered in the liter- 
ature (see Q for a review). A comprehensive list 
of possible sources with shock acceleration was 
thoroughly studied in ref.((T^]) with a conclusion, 
that maximum energy of acceleration does not ex- 
ceed 10 19 - f 20 eV (see also ref.(@) with a sim- 
ilar conclusion) . The most promising source from 
this list is a hot spot in radiogalaxy produced by 
a jet Jl£-21|, where maximum energy can reach 
~ 10 20 eV. Radiogalaxy M87, considered in p[], 
belongs to this class of sources. 

Gamma Ray Bursts (GRB) models offer two new 
mechanisms of acceleration to UHE. The first one 
p2| is acceleration by ultrarelativistic shock. A 
reflected particle gains at one reflection the en- 
ergy E ~ T 2 ,^, where T sh ~ 10 2 - 10 3 is the 
Lorentz factor of the shock and Ei is initial en- 
ergy of a particle. The second cycle of such ac- 
celeration has extremely low probability to occur 
^3|,Q and therefore to produce the particles with 
E ~ f 20 eV, this mechanism must operate in the 
space filled by pre-accelerated particles with en- 
ergies Ei > fO 14 eV. 



The second mechanism |25 works in the model 



with multiple shocks. The collisions of the shocks 
produces the turbulence where the particles are 
accelerated by Fermi II mechanism. The turbu- 
lent velocities are mildly relativistic in the fireball 
rest system. The maximum energy in the rest sys- 
tem , E' max ~ cHqI'q, is boosted by Lorentz factor 
r of fireball in laboratory system (here 1' and H' n 
are the maximum linear scale of turbulence with 
coherent magnetic field H' there). Taking for H' 
equipartition value, one obtains E ma . 



10 20 eV 



in the laboratory system. This mechanism faces 
two problems. Actually the maximum energy is 
somewhat less than 1 ~ 10 20 eV |pq| , if accel- 
eration time is evaluated more realistically. It 
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diminishes the energy of GZK cutoff in the dif- 
fuse spectrum, because it is formed by the parti- 
cles with production energies higher than the ob- 
served ones. The most serious problem, however, 
is that the produced flux of accelerated particles 
suffer the adiabatic energy losses Q . 

In summary, the acceleration (astrophysical) sce- 
narios are somewhat disfavoured, but not ex- 
cluded. Apart from them, many elementary par- 
ticle solutions were proposed to solve UHECR 
puzzle. Among them there is such an extreme 
proposal as breaking the Lorentz invariance P7[ , 
light gluino as the lightest supersymmetric parti- 
cle and UHE carrier |^8| , UHE neutrinos produc- 
ing UHECR due to resonance interaction with the 
dark matter neutrinos [ p9| and some other sugges- 
tions. In this paper I will review two most conser- 
vative sources of UHECR of non-accelerator ori- 
gin: Superheavy Dark Matter (SHDM) and Topo- 
logical Defects (TD). 

2. UHECR from Superheavy Dark Matter 

Superheavy Dark Matter (SHDM) as a source 
of UHECR was first suggested in refs.([^o|, 
|TJ). SHDM particles with masses larger than 
10 13 GeV are accumulated in the Galactic halo 
with overdensity - 10 5 and hence UHECR 
produced in their decays do not exhibit the GZK 
cutoff. The other observational signatures of this 
model are dominance of UHE photons ||(| and 
anisotropy connected with non-central position of 
the Sun in the Galactic halo ||^,|33|] . 

Production of SHDM 

SHDM particles are very efficiently produced 
by the various mechanisms at post-inflationary 
epochs. This common feature has a natural ex- 
planation. The SHDM particles due to their 
tremendous mass had never been in the thermal 
equilibrium in the Universe and never were rela- 
tivistic. Their mass density diminished as ~ 1/a 3 , 
while for all other particles it diminishes much 
faster as ~ 1/a 4 , where a is the scaling factor of 
the Universe. When normalized at inflationary 
epoch, aj = 1, a(t) reaches enormous value at 
large t. It is enough to produce negligible amount 
of superheavy particles in the post-inflationary 



epoch in order to provide fix ~ 1 now. Actually, 
in most cases one meets a problem of overproduc- 
tion of SHDM particles (further on we shall refer 
to them as to X-particles). 

One very general mechanism of X-particle pro- 
duction is given by creation of particles in time- 
variable classical field. In our case it can be 
inflaton field <f> or gravitational field. In case 
of inflaton field the direct coupling of X-particle 
(or some intermediate particle x) with inflaton is 
needed, e. g. g 2 <j> 2 X 2 or g 2 <fi 2 x 2 - The intermediate 
particle x then decays to X-particle. In case of 
time- variable gravitational field no coupling of X 
to inflaton or any other particles is needed: X- 
particles are produced due to their masses. For 
the review of above-mentioned mechanisms and 
references see Q. 

Super-heavy particles are very efficiently pro- 
duced at preheating p5[ . This stage, predeces- 
sor of reheating, is caused by oscillation of in- 
flaton field after inflation near the minimum of 
the potential. Such oscillating field can non- 
perturbatively (in the regime of broad paramet- 
ric resonance) produce the intermediate bosons 
X, which then decay to X-particles. The mass of 
X-particles can be one-two orders of magnitude 
larger than inflaton mass m^, which should be 
about 10 13 GeV to provide the amplitude of den- 
sity fluctuations observed by COBE. 

Another mechanism, more efficient than para- 
metric resonance and operating in its absence, is 
so-called instant preheating pq] . It works in the 
specific models, where mass of x particles is pro- 
portional to inflaton field, m x = g<p. When infla- 
ton goes through minimum of potential <fi = x~ 
particles are massless and they are very efficiently 
produced. When \(f>\ increases, m x increases too 
and can reach the value close to mp;. 

Another possible mechanisms of SHDM parti- 
cle production are non-equilibrium thermal pro- 
duction at reheating and by early topological de- 
fcct s^|30| . The latter can be produced at reheat- 

Gravitational production of particles occurs due 
to time variation of gravitational field during ex- 
pansion of the universe f37j| . For particles with 
the conformal coupling with gravity, (\/Q)RX, 
where R is the space-time curvature of the ex- 
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panding universe, the particle mass itself cou- 
ples a particle with the field (gravitation) and 
any other couplings are not needed. X particles 
can be even sterile! Neither inflation is needed 
for this production. It rather limits the gravi- 
tational production of the particles. Since this 
production is described by time variation of the 
Hubble constant H(t), only particles with masses 
mx < H (t) can be produced. In inflationary sce- 
nario H(t) < m0, where is the mass of the 
inflaton. It results in the limit on mass of pro- 
duced particles m x < 10 13 GeV 

The gravitational production of superheavy 
particles was recently studied in refs[ ^8|j39| (see 
p4| ] for a review). It is remarkable that for the 
mass mx ~ 10 13 GeV the relic density is fix ~ 1 
without any additional assumptions. It makes 
superheavy particles most natural candidates for 
Cold DM. 

Lifetime 

Superheavy particles are expected to be very 
short-lived. Even gravitational interaction (e.g. 
described by dimension 5 operators suppressed 
by the Planck mass) results in the lifetime much 
shorter than the age of the Universe to- The su- 
perheavy particles must be protected from fast 
decay by some symmetry, respected even by grav- 
itational interaction, and such symmetries are 
known. They are discrete gauge symmetries. 
They can be very weakly broken e.g. by wormholc 
effects J30| or instanton effects [[n) to provide the 
needed lifetime. The systematic analysis of bro- 
ken discrete gauge symmetries is given in ref . Eoj . 
For the group Ziq the lifetime of X-particle with 
mx ~ 10 13 — 10 14 GeV was found in the range 
10 11 — 10 26 yr. The realistic elementary parti- 
cle models for such long-lived particles were sug- 
gested |iH. 

Spectrum of UHECR 

Quark and gluons produced in the decay of su- 
perheavy particle originate QCD cascade, similar 
to that from Z° decay. The resulting spectrum 
of hadrons can be calculated using the standard 
QCD technique [ ^3p^ ]. The spectrum of hadrons 
is not power-law, its most spectacular feature is 
the Gaussian peak at small x. Photons dominate 
the primary spectrum by a factor ~ 6. The calcu- 



lations of the spectrum were performed in ref. |4q] 
(HERWIG MC simulation for ordinary QCD) and 
in ref.[||| (analytic MLLA calculations for SUSY 
QCD). 

Observational predictions. 

Overdensity 8 of SHDM particles in the Galactic 
halo is the same as for any other form of CDM, 
and numerically it is given by a ratio of CDM 
density observed in the halo to CDM density in 
extragalactic space (8 ~ 10 5 ). 

Spectrum of UHECR produced by decaying X- 
particles in the Galactic halo and beyond is shown 
in Fig. 2. One can see that UHE photon flux ap- 
preciably dominates over that of protons. 
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Figure 2. Predicted fluxes of UHE photons and 
protons from the decay of superheavy relic parti- 
cles with mass m x = 1 • 10 14 GeV. @/. The 
solid, dotted and dashed curves correspond to dif- 
ferent distributions of SHDM in the halo. Obser- 
vational data are from AG AS A. 

Anisotropy is caused by non-central position of 
the Sun in the halo. Most notable effect, the dif- 
ference in fluxes in directions of Galactic center 
and anticenter, cannot be observed by existing 
arrays. Calculated phase and amplitude of the 
first harmonic of anisotropy are compared 

in Fig. 3 with observations. In spite of the vi- 
sual agreement, one might only conclude that pre- 
dicted anisotropy does not contradict the obser- 
vations: within f.5cr AGASA data are compatible 
with isotropy. 
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estimate the fraction of UHE particles arriv- 
ing to us from one of these objects as / ~ 
(M c i/M h )(R 2 h /r 2 ), where R h ~ 100 kpc is a size 
of the halo, and M<j and Mh are the masses 
of a cluster and halo, respectively. For M c i ~ 
10 6 M©, M h - 10 12 M Q and r ~ 1 kpc, one ob- 
tains / ~ 0.01, i.e. about ten of such sources 
can provide the doublets and triplets observed 
in AGASA and other detectors. Part of these 
sources can be globular clusters. 

More detailed discussion will be presented in 
one of forthcoming publications by A.Vilenkin 
and the author. 



Figure 3. Amplitude and phase of the first har- 
monic of anisotropy for the AGASA (AG) and 
Yakutsk (YAK) arrays. Solid lines are for ISO 
distribution of DM and dots (NFW) for NFW nu- 
merical simulations J^H/ . BM and TW t 
calculations agree with each other. 



Angular clustering in UHECR arrival (doublet 
and triplet events) can be due to dumpiness of 
DM halo. Numerical N-body simulations show 
the presence of dense DM clouds in the halo. 
For example, the high resolution simulations of 
ref . |50| predict about 500 DM clouds with masses 
M ~ 10 8 M Q in the halo of our Galaxy. The 
baryonic content of these clouds should be low 
51 1, and therefore one cannot expect the iden- 
tification of all sources of UHECR doublets and 
triplets with the observed clouds. The smallest 
clumps resolved so far in the high resolution sim- 
ulations reach M c i ~ 10 6 M Q , i.e. they fall into 
range of the globular cluster masses. It could be 
that some of the doublet/triplet UHECR sources 
are globular clusters. 

The high resolution simulations demonstrate 
the early origin of the clumps (z w 5) (50), and 
therefore the core overdensity, as compared with 
present density, is (1 + z) 3 ~ 200. The extended 
halos of DM clouds can be stripped away by tidal 
interactions when clouds cross the galactic disc. 
The formation of dense compact DM objects was 
discussed in refs. [p2||5^] . 

Assuming the typical distance of a dense com- 
pact cloud to the sun as r ~ 1 kpc, one can 



3. Topological defects. 

Topological defects, TD, (for a review see [p4| ) 
can naturally produce particles of ultrahigh ener- 
gies (UHE). The pioneering observation of this 
possibility was made by Hill, Schramm and 
Walker []55[ (for a general analysis of TD as UHE 
CR sources see [pq] ). 

In many cases TD become unstable and de- 
compose to constituent fields, superheavy gauge 
and Higgs bosons (X-particles), which then de- 
cay producing UHECR. It could happen, for ex- 
ample, when two segments of ordinary string, or 
monopole and antimonopole touch each other, 
when electrical current in superconducting string 
reaches the critical value and in some other cases. 

In most cases the problem with UHECR from 
TD is not the maximal energy, but the fluxes. 
One very general reason for the low fluxes con- 
sists in the large distance between TD. A dimen- 
sion scale for this distance is the Hubble distance 
Hq 1 . However, in some rather exceptional cases 
this dimensional scale is multiplied to a small di- 
mensionless value r. If a distance between TD is 
larger than UHE proton attenuation length, then 
the flux at UHE is typically exponentially sup- 
pressed. 

The following TD have been discussed as 
potential sources of UHE particles: super- 
conducting strings |55|j , ordinary strings [f57f , 
]5q] , [j59| , magnetic monopoles , or more precisely 
bound monopolc-antimonopole pairs (monopolo- 
nium jS0|j6l] ] and monopole-antimonopole con- 
nected by strings networks of monopoles 
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connected by strings |6^] , necklaces (6^] , and vor- 
tons []65| . 

Monopolonia, monopole-antimonopole con- 
nected by strings and vortons are clustering in 
Galactic halo |Q and their observational signa- 
tures for UHECR are identical to SHDM particles 
discussed above. 

(i) Superconducting strings. 

As was first noted by Witten[|66|, in a wide 
class of elementary particle models, strings be- 
have like superconducting wires. Moving through 
cosmic magnetic fields, such strings develop elec- 
tric currents. Superconducting strings produce 
X particles when the electric current in the 
strings reaches the critical value. Superconduct- 
ing strings produce too small flux of UHE parti- 
cles |33| and thus they are disfavoured as sources 
of observed UHECR. 

(ii) Ordinary strings. 

There are several mechanisms by which ordinary 
strings can produce UHE particles. 

For a special choice of initial conditions, an or- 
dinary loop can collapse to a double line, releas- 
ing its total energy in the form of X-particles[^|. 
However, the probability of this mode of collapse 
is extremely small, and its contribution to the 
overall flux of UHE particles is negligible. 

String loops can also produce X-particles when 
they self-intersect (e.g. |67]]). Each intersection, 
however, gives only a few particles, and the cor- 
responding flux is very small (68| . 

Superheavy particles with large Lorentz fac- 
tors can be produced in the annihilation of cusps, 
when the two cusp segments overlap p9| ] . The en- 
ergy released in a single cusp event can be quite 
large, but again, the resulting flux of UHE parti- 
cles is too small to account for the observations 

It has been recently argued |)8| that long 
strings lose most of their energy not by produc- 
tion of closed loops, as it is generally believed, 
but by direct emission of heavy X-particles. If 
correct, this claim will change dramatically the 
standard picture of string evolution. It has been 
also suggested that the decay products of parti- 
cles produced in this way can explain the observed 
flux of UHECR [Q. However, as it is argued in 
ref |p3|, numerical simulations described in |58|| 



allow an alternative interpretation not connected 
with UHE particle production. 

But even if the conclusions of |5^| were correct, 
the particle production mechanism suggested in 
that paper cannot explain the observed flux of 
UHE particles. If particles are emitted directly 
from long strings, then the distance between UHE 
particle sources D is of the order of the Hubble 
distance Hq , D ~ Hq 1 3> R p , where R p is the 
proton attenuation length in the microwave back- 
ground radiation. In this case UHECR flux has 
an exponential cutoff at energy E ~ 3 ■ 10 10 GeV. 
In the case of accidental proximity of a string 
to the observer, the flux is strongly anisotropic. 
A fine-tuning in the position of the observer is 
needed to reconcile both requirements. 

(iii) Monopolonium and MM -pair connected by 
string. 

Monopole-antimonopole pairs (MM) can form 
bound state Q . Spiraling along the classical or- 
bits they fall to each other and annihilate, pro- 
ducing superheavy particles. The lifetime of this 
system depends on the initial (classical) radius, 
and can be larger than the age of the Universe to 
|i(J . Production of UHECR by monopolonia was 
studied in ref. |fil| (clustering of monopolonia in 
the Galactic halo was not noticed in this paper 
and was indicated later in ref. |}3| ) . 

Recently |62j it was demonstrated that fric- 
tion of monopoles in the cosmic plasma results in 
the monopolonium lifetime much shorter than to. 
Instead of monopolonium the authors have sug- 
gested a similar object, MM pair connected by 
a string, as a candidate for UHECR. This TD is 
produced in the sequence of the symmetry break- 
ing G — > H x U(l) — > H. At the first symmetry 
breaking monopoles are produced, at the second 
one each MM-pair is connected by a string. For 
the light strings the lifetime of this TD is larger 
than to- MM-pairs connected by strings are ac- 
cumulated in the halo as CDM and have the same 
observational signatures as SHDM particles. 

(iv) Network of monopoles connected by strings. 
The sequence of phase transitions 

G -> H x 1/(1) -> H x Z N (1) 

results in the formation of monopole-string net- 
works in which each monopole is attached to N 



7 



strings. Most of the monopoles and most of 
the strings belong to one infinite network. The 
evolution of networks is expected to be scale- 
invariant with a characteristic distance between 
monopoles d = nt, where t is the age of Universe 
and k = const. The production of UHE parti- 
cles are considered in |||. Each string attached 
to a monopole pulls it with a force equal to the 
string tension, u ~ rj 2 , where rj s is the symmetry 
breaking vev of strings. Then monopoles have a 
typical acceleration a ~ u/m, energy E ~ ad and 
Lorentz factor T m ~ fid/m, where m is the mass 
of the monopole. Monopole moving with accel- 
eration can, in principle, radiate gauge quanta, 
such as photons, gluons and weak gauge bosons, 
if the mass of gauge quantum (or the virtual- 
ly Q 2 in the case of gluon) is smaller than the 
monopole acceleration. The typical energy of ra- 
diated quanta in this case is e ~ Tmci- This en- 
ergy can be much higher than what is observed 
in UHECR. However, the produced flux (see |p3| ) 
is much smaller than the observed one. 
(v) Vortons. 

Vortons are charge and current carrying loops of 
superconducting string stabilized by their angu- 
lar momentum JnJ. Although classically stable, 
vortons decay by gradually losing charge carriers 
through quantum tunneling. Their lifetime, how- 
ever, can be greater than the present age of the 
universe, in which case the escaping A-particles 
will produce a flux of cosmic rays. The X-particle 
mass is set by the symmetry breaking scale rjx of 
string superconductivity. 

The number density of vortons formed in the 
early universe is rather uncertain. According to 
the analysis in ref.j72|, vortons are overproduced 
in models with r/x > 10 9 GeV, so all such mod- 
els have to be ruled out. In that case, vortons 
cannot contribute to the flux of UHECR. How- 
ever, an alternative analysis Q] suggests that 
the excluded range is 10 9 GeU < n x < !0 12 GeV, 
while for r\x ^> 10 12 GeV vorton formation is 
strongly suppressed. This allows a window for 
potentially interesting vorton densities with rjx ~ 
10 12 - 10 13 GeV. Production of Ultra High En- 
ergy particles by decaying vortons was studied in 
ref.@. 

Like monopoles connected by strings and SH 



relic particles, vortons are clustering in the Galac- 
tic halo and UHECR production and spectra are 
identical in these three cases. 

(vi) Necklaces. 
Necklaces are hybrid TD corresponding to the 
case N = 2 in Eq. ([!]), i.e. to the case when each 
monopole is attached to two strings. This system 
resembles "ordinary" cosmic strings, except the 
strings look like necklaces with monopoles play- 
ing the role of beads. The evolution of necklaces 
depends strongly on the parameter 



r = mj [id, 



(2) 



where d is the average separation between 
monopoles and antimonopolcs along the strings. 
As it is argued in ref . , necklaces might evolve 
to configurations with r 1, though numeri- 
cal simulations are needed to confirm this con- 
clusion. Monopoles and antimonopoles trapped 
in the necklaces inevitably annihilate in the end, 
producing first the heavy Higgs and gauge bosons 
(X-particles) and then hadrons. The rate of X- 
particle production can be estimated as 64 
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t 3 mx 



(3) 



Restriction due to e-m cascade radiation de- 
mands the cascade energy density lu ccis < 2 ■ 
10~ 6 eV/cm 3 . The cascade energy density pro- 
duced by necklaces can be calculated as 

rta dt 1 _ 3 2 fj, 



1 



" ' ■ < " . J - T' I* 

1 JO 



t 3 (1 



(4) 



(i 



where f v ~ 0.5 is a fraction of total energy re- 
lease transferred to the cascade. The separation 
between necklaces is given by |34| D ~ r -1 / 2 ^ 
for large r. Since r 2 u is limited by cascade radi- 
ation, Eq.(^), one can obtain a lower limit on the 
separation D between necklaces as 

1/4 

D - [ .T' lf ~ 1 t > 10(/V10 6 GeV 2 ) 1/4 kpc, (5) 



Thus, necklaces can give a realistic example of 
the case when separation between sources is small 
and the Universe can be assumed uniformly filled 
by the sources. 

The fluxes of UHE protons and photons are 
shown in Fig. 4 according to calculations of 
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ref. |Q. Due to absorption of UHE photons 
the proton-induced EAS from necklaces strongly 
dominate over those induced by photons at all 
energies except E > 3 • 10 11 GeV, where photon- 
induced showers can comprise an appreciable 
fraction of the total rate. 
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Figure 4. Spectrum of UHE photons + protons 
from necklaces compared with AG AS A measure- 
ments. Solid and dashed curves correspond to dif- 
ferent absorption of UHE photons in extragalactic 
space. 

4. Conclusions 

At E > 1 • 10 19 eV a new component of cosmic 
rays with a flat spectrum is observed. According 
to the Fly's Eye and Yakutsk data the chemical 
composition is better described by protons than 
heavy nuclei. The AG AS A data are consistent 
with isotropy in arrival of the particles, but about 
20% of particles at E > 4 • 10 19 eV arrive as 
doublets and triplets within ~ 2 — 4° . 

The galactic origin of UHECR due to conven- 
tional sources is disfavoured: the maximal ob- 
served energies are higher than that calculated 
for the galactic sources, and the strong Galactic 
disc anisotropy is predicted even for the extreme 
magnetic fields in the disc and halo. 

The signature of extragalactic UHECR is GZK 
cutoff. The position of steepening is model- 
dependent value. For the Universe uniformly 
filled with sources, the steepening starts at Ebb ~ 
3 • 10 19 eV and has E 1/2 « 6 • 10 19 eV (the energy 
at which spectrum becomes a factor of two lower 



than a power-law extrapolation from lower ener- 
gies) . The spectra of UHE nuclei exhibit steepen- 
ing approximately at the same energy as protons. 
UHE photons have small absorption length due 
to interaction with radio background radiation. 

The extragalactic astrophysical sources theo- 
retically studied so far, have either too small 
E m ax or are located too far away. The Local 
Supcrcluster (LS) model can give spectrum with 
E 1 / 2 ~ 10 20 eV, if overdensity of the sources is 
larger than 10. However, IRAS galaxy counts give 
overdensity S = 1.4. 

GRBs and a nearby single source (e.g. M87) 
remain the potential candidates for the observed 
UHECR. 

Superheavy Dark Matter can be the source 
of observed UHECR. These objects can be relic 
superheavy particles or topological defects such 
as MM-pairs connected by strings or vortons. 
These objects are accumulated in the halo and 
thus the resulting spectrum of UHECR does not 
have the GZK cutoff. In this case UHECR is 
a signal from inflationary epoch, because both 
superheavy particles and topological defects are 
most probably produced during reheating. 

The observational signatures of UHECR from 
SHDM are (i) absence of GZK cutoff, (ii) UHE 
photons as the primaries and (iii) anisotropy due 
to non-central position of the Sun in the halo. 
The angular clustering is possible due to dumpi- 
ness of DM in the halo. 

Topological Defects naturally produce particles 
with extremely high energies, much in excess of 
what is presently observed. However, the fluxes 
from most known TD are too small. Only neck- 
laces, MM connected by strings and vortons re- 
main candidates for the sources of the observed 
UHECR. Necklaces give so far the only known 
example of extragalactic TD as a sources of 
UHECR. Its signature is the presence of the pho- 
ton component in the primary radiation and its 
dominance at the highest energies E > 10 20 eV. 
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